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Abstract 

In the framework of analytic approach to QCD, which has been recently intensively developed, 
the dependence of nonperturbative contributions in a running couphng of strong interaction on ini- 
tial perturbative approximation to 3-loop order is studied. The nonperturbative contributions are 
obtained in an explicit form. In the ultraviolet region they are shown to be represented in the form of 
the expansion in the inverse powers of Euclidean momentum squared. The expansion coefficients are 
calculated for different numbers of active quark flavors n/ and for different numbers of loops taken 
into account. For aU n/ of interest it is shown that 2-loop order and 3-loop order corrections result in 
partial compensation of 1-loop order leading in the ultraviolet region nonperturbative contribution. 
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It is widely believed that unphysical singularities of the perturbation theory in the infrared region 
of QCD should be canceled by the nonperturbative contributions. The nonperturbative contributions 
arise quite naturally in an analytic approach [Q to QCD. The so-called " analyticization procedure" is 
used in this approach. The main purpose of this procedure is to remove nonphysical singularities from 
approximate (perturbative) expressions for the Green functions of QFT. The idea of the procedure goes 
back to Refs. ||] devoted to the ghost pole problem in QED. The foundation of the procedure is 
the principle of summation of imaginary parts of the perturbation theory terms. Then, the Kallen - 
Lehmann spectral representation results in the expressions without nonphysical singularities. In recent 
papers || it is suggested to solve the ghost pole problem in QCD demanding the as{q^) be analytical 
in (to compare with dispersive approach ||^). As a result, instead of the one-loop expression ai^-* (g^) ~ 
(47r/6o)/ ln((j'^/A^) taking into account the leading logarithms and having the ghost pole at = {q"^ 
is the Euclidean momentum squared), one obtains the expression 
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In(q2/A^) A2-g: 
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Eq. (|l|) is an analytic function in the complex g^-plane with a cut along the negative real semiaxis. The 
pole of the perturbative running coupling at q^ = is canceled by the nonperturbative contribution 
(A^ |g2_^Q ~ /i^ exp{— (47r)2/(6o(;2)|) g^^^d the value ain(O) = 47r/6o appeared finite and independent 
of A. The most important feature of the "analyticization procedure" discovered |^ is the stability 
property of the value of the " analytically improved" running coupling constant at zero with respect to 
high corrections, Q;in(0) = Q;in(0) = Q;in(0). This property provides the high corrections stabihty of 
OLani'f') in the whole infrared region. 

The 1-loop order nonperturbative contribution in Eq. (|l|) can be presented as convergent at q^ > A^ 
constant signs series in the inverse powers of the momentum squared. For "standard" as well as for the 
iterative 2-loop perturbative input the nonperturbative contributions in the analytic running coupling are 
calculated explicitly in Ref. In the ultraviolet region the nonperturbative contributions can be also 
represented as a series in inverse powers of momentum squared. In this paper we extract in an explicit form 
the nonperturbative contributions to aanil'^) up to the 3-loop order in the analytic approach to QCD, 
and investigate their ultraviolet behavior. To handle the singularities originating from the perturbative 
input, we develop here the method which is more general than that of Ref. 0. 

According to the definition the analytic running coupling is obtained by the integral representation 

oo 

1 f da 

aan(x) = - / ■ p(cr), (2) 

n J X -\- a 



where the spectral density p{a) = Imaa„(— cr — iO). It is seen that dispersively-modified coupling of 
form has an analytical structure which is consistent with causality. According to the analytic approach 
to QCD we adopt that ImaQ„(— ct — iO) = Im a(—o' — iO) where a(x) is an appropriately normalized 
perturbative running coupling. 

The behavior of the QCD running coupling as{fi^) is defined by the renormalization group equation 

fi^^4^ ^ P{a,) ^ Po^l + Pial + P^c^i + (3) 



where the coefficients 



1 2 

Po = -^^^0, foo = 11 - -^^nf, 

ZTT O 

1 19 

/^i = -^^i' bi=51-—nf, 

1 . . nor^ 5033 325 2 

02 = -^62, b2 = 2857 - -^n, + —n}. (4) 

The ffi'st two coefficients /?o, /3i do not depend on the renormalization scheme choice. Here Uf is the 
number of active quark flavors. The standard three- loop solution of Eq. (j^) is written in the form of 
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expansion in inverse powers of logarithms 



bo ln(/i2/A2 



1 - 



2bi ln[ln(/iVA2)] 
1^ ln(AiVA2) 
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Let us introduce the function a{x) — {bo/4n)as{q'^), where x = g^/A^. Then, instead of one can write 

'ln2(lnx) ln(lna;) n 



1 hi Inx) 

a{x) = b—^ 

m X In"' X 

where the coefficients b and k are equal to 
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At Uf = 3 bg = 9, and b — 64/81 ~ 0.7901, k ~ 0.4147. At x ~ 1 the perturbative running coupling is 
singular. At large x the 1-loop term of Eq. (||) defines the ultraviolet behavior of a{x) but for small x 
the behavior of the running coupling depends on the approximation we adopt, and at a; = 1 there are 
singularities of a different analytical structure. Namely, at a; ~ 1 



a'^^Hx) 
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This is not an obstacle for the analytic approach which removes this nonphysical singularities. By making 
the analytic continuation of Eq. (o) into the Minkowski space x = —a — iO, one obtains 



a{—a — iO) = 
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in (In (7 — i7r)2 
ln(ln a — in) 
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Function a{x) in Eq. (|^) is regular and real for real x > 1. So, to find the spectral density p{a) we can 
use the reflection principle {a{x))* = a{x*) where x is considered as a complex variable. Then 
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(a(-CT - iO) - a{~(7 + iO)) . 



By the change of variable of the form a — exp(i), the analytical expression is derived from 
as follows: 
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[t - in)^ {t + in)^ (t - in)'-^ [t + in)^ _ 

Let us see what the singularities of the integrand of (^l]) in the complex t-plane are. First of all the 
integrand has simple poles at i = Ina; ± in{l + 2n), n — 0,1,2,.... All the residues of the function 
exp(t)/(a; + exp(t)) at these points are equal to unity. Apart from these poles the integrand of (|ll|) has 
simple poles at i = ±i7r, the third order poles and logarithmic type branch points which coincide with 
the second order and third order poles. Let us cut the complex i-plane in a standard way, t = ±«7r — A, 
with A being the real parameter varying from to oo. The integrand in (fO) multiplied by t goes to 
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zero at I i |— > oo. That allows one to append the integration by the arch of the "infinite" radius without 
affecting the value of the integral. Close the integration contour Ci in the upper half-plane of the complex 
variable t excluding the singularity at t = in. In this case an additional contribution emerges due to the 
integration along the sides of the cut and around the singularities at t = in. The corresponding contour 
we denote as C2. 

Let us turn to the integration along contour Ci . For the integrand of Eq. (|l^) which we denote as 
F{t) the residues at t = lii{x) + in{l + 2n), n = 0, 1, 2, ... are as follows 
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(ln(a;) + 2nin)^ {\n{x) + 2m{n + l))^ (ln(a:) + 27rm)3 (ln(a;) + 2m{n + 1))^ 

By using the residue theorem one readily obtains the contribution A(x) to integral ( pj]) from the inte- 
gration along contour Ci. It reads 
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F{t) dt 
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Res F{t = ln(x) + i7r(l + 2n)) 
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We can see that this contribution is exactly equal to initial Eq. (^. Therefore, we call it a perturbative 
part of Qanix), aP*{x) = A(a;). The remaining contribution of the integral along contour C2 can be 
naturally called a nonperturbative part of aan (x) , 



aan{x)^aP\x)+aZ'{x). 

Let us turn to the calculation of a^JJ*(a:). We can omit the terms of the integrand in Eq. 
no singularities at t = in. Then, we have 
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ll) which have 



Let us change the variable t 



\n'{t — in) \n{t — in) ^ k 
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in and introduce the function 
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1 — X exp(— z) 



Then, we can rewrite Eq. ( |15D in the form 
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The cut in the complex z-plane goes now from zero to — cx3. Starting from z = —00 — iO contour C goes 
along the lower side of the cut, then goes around the origin, and next it goes further along the upper side 
of the cut to z = —00 + iO. Here we consider x as a real variable, x > 1. Then, contour C can be chosen 
in such a way that it does not envelop "superfluous" singularities and the conditions used in Appendix 
to find the corresponding integrals be satisfied. Function (n™ with its derivatives decrease exponentially 
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at z — cx), therefore, we shall omit the boundary terms in formulas given in Appendix. We shall need 
further the explicit expressions for the derivatives of f{z), which read 



X exp(— z) 



/'"(^) = - 



(1 — X exp(— ' 
X cxp(— z) 



X cxp(— z)(l + X cxp(— z)) 



(1 — a;exp(— z))' 



(1 — a; exp(— z))"^ 
(l + Ax exp(-z) + x'^ exp(-2z)) 



From Eq. (0) one can obtain 

<ni^) = -^^ I dz |/'(z)ln(z)-6/"(z) {\n{z) + \\n\z)^ + 



(18) 



1 + -A^)ln(z) + -In2(z) + -In3(z) 
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Taking into account that function /(z) with its derivatives is regular at real negative semiaxis of z, we 
can rewrite equation (|l9[) in the form 



du 



f'iu)A,{u) - bf"{u) [ A^{u) + ^A,iu) ] + 



+ b^riu) + -Kj A^{u) + -A^{u) + -A3(u 
where u is real, m < and Aj(w) are discontinuities of the logarithms 
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Ai{u) = — (ln(u + iO)-ln(M-iO)) = 1, 

A2{u) = ^ (ln^(M + iO) -ln^(M-iO)) = 21n(-u), 

A3(w) = — (ln^(u + iO) -ln^(u-iO)) = 31n^(-u) -TT^ 
27ri 



Let us introduce the variable a — exp(u). From Eqs. (|18|), (poD, ( |21| ) we obtain 
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aZ\x) = -x / da 



{x — (t)2 



6[l + ln(-ln(a))] 



X + a 



{x — a)'^ 
x^ + Axa + cr^ 



{x - a)4 

Integrating the terms of Eq. (^) independent of logarithms one can obtain 
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K TT^ \ a;(l + x) 
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da 



ln(-ln(a)) + iln2(-ln(^)) 



x^ + Axa + a^ 

(2; - (7)4 



(23) 



This formula gives the nonperturbativc contributions in an explicit form and is convenient for numerical 
study. 



4 



Let us turn to the ultraviolet behavior of the nonperturbative contributions. Expanding Eq. ( [23| ) in 
inverse powers of x we have 

oo 

a/' 



npt 



(24) 



where 



-1 + bn 



1 + ny dcrcr"-Mn(-ln(cr)) 








ln(-ln(a)) + iln2(-ln(a)) 
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Making the change of variable a = exp(— i) and integrating ||T^ over t, one can find 



da cr"-i In (- In(cr)) = / dt e""* hi(t) = (ln(n) + 7) , 

J n 




(26) 



d(7(7"-Mn'(-ln((T)) = 

From Eqs. M), M), (|^ we finally have 



(ln(n) + 7)' + ^ 
6 



(27) 



-1 + 6n (1 ~ 7 - ln(n)) - 



6 
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Here 7 is the Euler constant, 7 ~ 0.5772. We can see from Eq. ( p8[ ) that power series ( p4|) is uniformly 
convergent at a; > 1 and its convergence radius is equal to unity. To see the details of the coefficients 
Cn behavior we give the results of the numerical study in Table I. The coefficients c„ depend on n, n/ 
and the number of loops taken into account. The 1-loop order contribution to Cn equals —1 for all n and 
Uf. The 2-loop order corrections to c„ decreases monotonously with increasing n for Uf = 0,3,4,5,6. 
In this case the correction to the leading term (n = 1) is positive (it varies from 0.36 at n/ = to 0.22 
at n/ = 6), whereas for the next terms the corrections are negative. The 3- loop order corrections to c„ 
increases at first and then monotonously decreases with increasing n for all Uf we consider. 

In the ultraviolet region (cc ^ 1) the nonperturbative contributions are determined by the first term 
of the series (^4|). At Uf < 4 the 3-loop order corrections change insignificantly the 2-loop order results 
for nonperturbative contributions. At n/ = 5 the 3-loop order correction is 2.4 times less than the 2-loop 
order correction and at n/ = 6 it exceeds the 2-loop order correction. We also note that with increasing 
rif the decrease of the 2-loop corrections is equilibrated to some excess by the increase of the 3-loop 
corrections. 

It seems to be most interesting that an account for the 2-loop order corrections results in some 
compensation of the 1-loop order leading at large x term of the form 1 /x and at 3-loop order compensation 
remains valid. We can see from Table I that at Uf > 3 the 3-loop order |ci| is somewhat less than 2-loop 
order |ci|. However, it apparently does not mean that in an analytic approach the high-loop corrections 
lead to the total compensation of nonperturbative contributions in the ultraviolet region. 

I am indebted to B.A. Arbuzov, V.A. Petrov, V.E. Rochev for useful discussions. This work has been 
supported in part by RFBR under Grants No. 99-01-00091, No. 98-02-16690. 

Appendix 

We give here the identities we need in our computations which can be obtained by means of an 
integration by parts. Let function f{z) of complex variable z is regular in some domain D where z — 
& D. Dealing with the singularities of the integrands at the origin of the pole type coinciding with the 
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logarithmic type branch points, we cut domain D along real negative semiaxis. Then, for any contour C 
in cut domain D, which goes from zi ^ to Z2 0, one can find 



dz\n{z)f{z)+\n{z)f{z) , 

Zl 

J $-^(^) = ~^ / dzHz)f"{z) + i^-^f{z)-^J'{z) + ^lniz)f"{z)^ 

c c 

I J ln(^)/(^) = -ldz (ln{z) + ^ ]n\z)^ f"{z) + [-\f{z) - ^/(^) + ln(^)/'(^) + 



+ 2!^ (^)/ 



ln(z) 



j ^Hz)m = -jdz gln(z) + \\n\z)^ f"{z) + {-^/W - ^/{z) - '-^f{z) - 



2^2 



/ pAz)m = -l dz (JlnW + ^ln^(.) + iln^(.))r(.)+{-^/(.)-^/'(.)-^/(.)- 



2z 



2z 
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Table 1: The dependence of c„ on n and n/ for 1-loop, 2-loop and 3-loop cases. 





n 


^\—loop 
n 


^2-loop 


^3 -loop 


J2—loop 
n 


^3— loop 
n 


n/ = 


1 


-1.0 


0.35640 


-0.01568 


-0.64360 


-0.65929 




2 


-1.0 


-0.45582 


0.08741 


-1.45582 


-1.36841 




3 


-1.0 


-1.70912 


-1.03012 


-2.70912 


-3.73924 




4 


-1.0 


-3.24886 


-4.51236 


-4.24886 


-8.76122 




5 


-1.0 


-5.00160 


-11.31238 


-6.00160 


-17.31398 




6 


-1.0 


-6.92407 


-22.24971 


-7.92407 


-30.17378 




7 


-1.0 


-8.98770 


-38.04599 


-9.98770 


-48.03369 




8 


-1.0 


-11.17217 


-59.34790 


-12.17217 


-71.52007 




9 


-1.0 


-13.46228 


-86.74274 


-14.46228 


-101.20502 




10 


-1.0 


-15.84626 


-120.76945 


-16.84626 


-137.61571 


n/ = 3 


1 


-1.0 


0.33405 


0.01608 


-0.66595 


-0.64987 




2 


-1.0 


-0.42724 


0.19623 


-1.42724 


-1.23101 




3 


-1.0 


-1.60196 


-0.63626 


-2.60196 


-3.23823 




4 


-1.0 


-3.04517 


-3.48652 


-4.04517 


-7.53168 




5 


-1.0 


-4.68801 


-9.19185 


-5.68801 


-14.87987 




6 


-1.0 


-6.48996 


-18.47225 


-7.48996 


-25.96221 




7 


-1.0 


-8.42420 


-31.96169 


-9.42420 


-41.38590 




8 


-1.0 


-10.47171 


-50.22832 


-11.47171 


-61.70003 




9 


-1.0 


-12.61824 


-73.78810 


-13.61824 


-87.40634 




10 


-1.0 


-14.85275 


-103.11452 


-15.85275 


-118.96726 


n/ = 4 


1 


-1.0 


0.31252 


0.04949 


-0.68748 


-0.63799 


2 


-1.0 


-0.39970 


0.31340 


-1.39970 


-1.08630 




3 


-1.0 


-1.49872 


-0.23818 


-2.49872 


-2.73690 




4 


-1.0 


-2.84891 


-2.48499 


-3.84891 


-6.33389 




5 


-1.0 


-4.38587 


-7.15989 


-5.38587 


-12.54576 




6 


-1.0 


-6.07168 


-14.89306 


-7.07168 


-21.96474 




7 


-1.0 


-7.88126 


-26.23939 


-8.88126 


-35.12065 




8 


-1.0 


-9.79681 


-41.69613 


-10.79681 


-52.49294 




9 


-1.0 


-11.80500 


-61.71490 


-12.80500 


-74.51990 




10 


-1.0 


-13.89549 


-86.71011 


-14.89549 


-101.60561 


n/ = 5 


1 


-1.0 


0.27813 


0.11653 


-0.72187 


-0.60535 


2 


-1.0 


-0.35571 


0.55755 


-1.35571 


-0.79817 




3 


-1.0 


-1.33377 


0.50736 


-2.33377 


-1.82641 




4 


-1.0 


-2.53536 


-0.73077 


-3.53536 


-4.26613 




5 


-1.0 


-3.90317 


-3.73728 


-4.90317 


-8.64045 




6 


-1.0 


-5.40344 


-9.01126 


-6.40344 


-15.41470 




7 


-1.0 


-7.01386 


-16.99218 


-8.01386 


-25.00604 




8 


-1.0 


-8.71859 


-28.07387 


-9.71859 


-37.79246 




9 


-1.0 


-10.50576 


-42.61401 


-11.50576 


-54.11977 




10 


-1.0 


-12.36618 


-60.94080 


-13.36618 


-74.30698 


n/ = 6 


1 


-1.0 


0.22433 


0.25378 


-0.77567 


-0.52189 




2 


-1.0 


-0.28692 


1.07460 


-1.28692 


-0.21231 




3 


-1.0 


-1.07581 


1.93179 


-2.07581 


-0.14402 




4 


-1.0 


-2.04500 


2.37205 


-3.04500 


-0.67296 




5 


-1.0 


-3.14826 


2.01775 


-4.14826 


-2.13052 




6 


-1.0 


-4.35837 


0.54419 


-5.35837 


-4.81418 




7 


-1.0 


-5.65732 


-2.33454 


-6.65732 


-8.99186 




8 


-1.0 


-7.03234 


-6.87466 


-8.03234 


-14.90700 




9 


-1.0 


-8.47386 


-13.30888 


-9.47386 


-22.78274 




10 


-1.0 


-9.97445 


-21.85073 


-10.97445 


-32.82518 
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